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ABSTRACT 

We describe a new method that only focuses on the local region surrounding an infalling dwarf in an effort to understand how the 
hot baryonic halo alters the chemodynamical evolution of dwarf galaxies. Using this method, we examine how a dwarf, similar to 
Sextans dwarf spheroidal, evolves in the corona of a Milky Way-like galaxy. We find that even at high perigalacticons the synergistic 
interaction between ram pressure and tidal forces transform a dwarf into a stream, suggesting that Sextans was much more massive 
in the past to have survived its perigalacticon passage. In addition, the large confining pressure of the hot corona allows gas that was 
originally at the outskirts to begin forming stars, initially forming stars of low metallicity compared to the dwarf evolved in isolation. 
This increase in star formation eventually allows a dwarf galaxy to form more metal rich stars compared to a dwarf in isolation, 
but only if the dwarf retains gas for a sufficiently long period of time. In addition, dwarfs that formed substantial numbers of stars 
post-infall have a slightly elevated [Mg/Fe] at high metallicity ([Fe/H]—1.5). 
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1. Introduction 

The concordant cosmological model of dark energy and cold 
dark matter, ACDM, predicts that galactic halos contain numer- 
ous subhalos that have accreted over cosmic time (e.g. |Klypin| 
et al.| 1999} [Moore et al.|1999| ). A portion of these accreted sub¬ 
halos are thought to host the dwarf galaxies that are visible today 
as satellite galaxies, although exactly which subhalos would host 
galaxies is an open question fBoylan-Kolchin et al.|2011| ). 

Regardless of the details of early galaxy formation, only the 
most massive satellite dwarf galaxies orbiting close in to the host 
galaxy are able to retain gas to the present day ( [Grcevich & Put-| 
man|2009t [Chiboucas et al.||2013| ), with the other dwarfs being 
gas-poor dwarf spheroidals (dSphs). Further out, however, dwarf 
galaxies of lower mass are able to retain a detectable amount of 
gas, which is able to fuel ongoing star formation, suggesting that 
environmental effects are the cause of this gas-poor gas-rich di¬ 
chotomy ( [Grcevich & Putman|2009| ). 

As galaxy simulations increase in resolution more attention 
is being paid towards the evolution of dwarf galaxies and how 
simulations compare to observations. Previously we have ex¬ 
amined the chemodynamical ev olution of dSphs as t hey inter¬ 
act tidally with a larger galaxy ( [Nichols et al.||2014l hereafter 
NRJ14). Considering only tidal effects, dwarfs lose substantial 
quantities of gas because of a synergistic effect of supernova in¬ 
duced expulsion and tidal forces. However, all dwarfs were able 
to retain quantities of gas that would be readily detectable today 
unless they were completely destroyed by tides. A potential solu¬ 
tion to the removal of gas, without tidally destroying the dwarfs, 
lies in ram pressure stripping. 

Previous work has shown that the interaction of the hot 
corona and satellite galaxies is likely to remove all gas from them 


[et al.|2006l Zolotov et al.|2012[), or the treatment of star forma- 


however, problems with Kelvin-Helmholtz instabilities ( [Mayer| 


tion (Nichols & Bland-Hawthorn|2Qll][2013||Gatto et al.|2Q13[ ) 

justify new simulations. 

The stripping of galaxies close in (< 50 kpc) is likely to 
occur independent of any internal process within the galaxy 
( [Mayer et al.|2006[ ), but further out, internal feedback resulting 
from star formation is needed to remove the gas ( [Gatto etaL] 
[2013[ ). In previous work on the topic, the star formation has only 
been treated with simplistic models for analytic ease ( [Nichols ^ 
[Bland-Hawthom|201 1[[2013] ) or assu med based on the o bserva- 
tionally derived star formation history ( [Gatto et al.|2013[ ). These 
approaches gloss over the wide range of observational properties 
that exist within the dSphs and what that might say about their 
evolution post-infall. 


[Sawala et al.[ ( [2012[ ) and [Zolotov et al.[ ( [2012[ ) simulated 
dwarf galaxies within a cosmological framework, running sim¬ 
ulations of Milky Way sized galaxies at high resolution. Here, 
using traditional smoothed particle hydrodynamic (SPH) simu- 
latio ns with self-co nsistent star formation and metallicity injec- 
tion, [Zolotov et al. [ found that baryons can assist in the tidal strip¬ 
ping of dwarf galaxies and can dramatically change the inner re- 
gion of a dwarf through supernova feedb ack (see also [Governato 
et al.|2010{ [Pontzen & Governato|2012|). Kno wn problems with 


the traditional SPH models ( [Agertz et al.[2007[ ) however, allowed 
many dwarf spheroidal analogues to retain gas to the present day 
and potentially underestimated the impact of ram pressure strip¬ 
ping. 


([Mayer et al.[[2006[ [N ichols & Bl a nd-Hawth or nj[2011t Zolotov[ 

[et al.|2012|[Nichols & Bland-Hawthom[2013[]Gatto et al. 2013[), 


* e-mail: matthew. nichols@epfl. ch 


[Sawala et al.[ found that many dwarfs are likely to enter the 
halo already gas-free, and tidal stripping combined with UV 
heating and supernova feedback can remove gas from many of 
the remaining dwarfs. With the use of traditional SPH models 
and at low resolution, however, ram pressure stripping may have 
been overestimated for some dwarfs, which may accelerate the 
stripping and cause some dwarfs to be tidally destroyed through 
this baryon-tide synergy. 
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Ideally large cosmological simulatio ns, such as found in 
|Sawala et al.| ( 2Q12| ); |ZolotW et al. (20121, could be run for any 
change in the initial conditions or numerical methods (e.g. mesh¬ 
less or updated SPH methods). Because of the large computa¬ 
tional cost, however, more efficient methods must be considered 
for the bulk of simulations. 

The need for star formation, metallicity evolution, and tidal 
forces to be handled simultaneously can be understood by con¬ 
sidering the wide range of properties impacted by internal star 
formation and tidal or ram-pressure stripping. For exam ple, 
present day gas-rich dwarfs show constant density cores ( |Oh | 
et al.|2015]), a process thought to be caused by superno va feed- 


back ( [Governato et al.|2010t[Pontzen & Govemato|2012| ). On the 
scales of the dwarf spheroidals the existence of cores is slightly 
less clear (see e.g. |de Blok et al.||2008[ [Walker & Penarrubla] 
2011[[Breddels & Helmi|2Q13|) and, as the inner profile can alter 


the tidal evolution ( [Kazantzidis et al.|2013| ), simulating the inter¬ 
nal star formation is vitally important to understanding the evo¬ 
lution of a dwarf. In addition to the internal star formation chang¬ 
ing its profile, the removal of gas also impacts the inner density 
( Arrakiet al.|2014| ) and the inclusion of this baryon loss may lead 
to a lower fraction of surviving dwarfs than expected from sim¬ 
ulations that retain gas. In addition the metallicity evolution is 
intricately linked with the star formation history (and contributes 
towards its regulation through cooling), and with high-resolution 
studies of dwarf spheroidals now available (e .g. [Shetrone et slI 


lOOTtlAoki et al.|2009tfiMelmeyer et al.|20T0l[Kirby et al 2011 


Theler et al. |2Q 15 [ Just for Sextans dwarf spheroidal) the metal¬ 

licity evolution can provide a sanity test for any model aiming to 
reproduce small satellite galaxies. 

In this paper, we describe modifications to the chemody- 
namical tree-SPH code GEAR ( [Revaz & Jablonka[[20I2| ), which 
allows us to efficiently simulate the post-accretion history of a 
galaxy and discuss how a dwarf galaxy, which has properties 
similar to an observed dwarf (excepting of course the excess in 
gas) changes as it evolves into a host galaxy after infall. We also 
pay particular attention to the chemical evolution of the dwarf 
galaxies and how this evolution is influenced by ram pressure 
stripping. We describe the method in ^ the initial conditions 
and set up of our simulations in ^ the post-infall evolution by 
comparison with an isolated evolution in ^ and conclude in ^ 


2. Simulating ram pressure stripping 

We have previously used the chemodynamical code GEAR iEI 
vaz & Jablonka[[2012|), which is an extension of the public 
Gadget-2 code ^Springel|20Q5] ), to examine dwarfs in isolation 
IRevaz & Jablonka[[2012| ) andthose that are tidally interacting 
|NRJ14[ ). Along with the latest additions to the code ( [Revaz et aT 


2Q15|), including the pressure-entropy SPH formulation of Hop 


kins[ ( [201^ , we extend this code to simulate the ram pressure 
and tidal forces simultaneously upon an infalling dwarf. 

To avoid the computationally expensive simulation of an en¬ 
tire hot halo, we performed wind-tunnel simulations, where only 
the region of the halo that interacts with the dwarf is considered. 
The dwarf is placed inside a gas-filled periodic box, where par¬ 
ticles are injected from one side and deleted when crossing the 
opposite side. The use of the pressure-entropy SPH for mulation 
allows us to account for Kelvin-Helmholtz instabilities ([Hopkins 


2013|) a kno wn problem with the traditional SPH models ^gertz 


et al.|2007| ). To properly describe the tidal forces acting upon the 
dwarf, we rotate this periodic box along the orbit of the dwarf. 
This kind of approach allows us to capture any synergy between 
the tidal and ram pressure forces upon the dwarf. The moving 
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of the periodic box is conc eptually similar to the method used 
by [Ploeckinger et aT ( [2014 ) to study the evolution of tidal dwarf 
galaxies, but with an evolving halo density. We describe this ro¬ 
tation scheme in §2.1[ 

In addition to the tidal forces changing over the orbit of the 
dwarf, the changing density of the galactic corona should be re¬ 
flected in the density of gas within the periodic box. In addition, 
the wake behind a dwarf should not be able to cross the periodic 
boundary to prevent it from interacting with the dwarf. These 
issues can be addressed through the deletion of particles at the 
boundary behind the dwarf, and an insertion of new gas in front 
of the dwarf. This kind of approach has already been used suc- 
cessfu lly in grid codes to study the hot halo density ( [Gatto et ^ 


[201 3[) and the wake prod uced by galaxies infalling into a clus- 
ter ( [Roediger et al.|2015[ ). We describe a similar adaption for the 
SPH code GEAR here, described in |2.2[ and §2.3[ for the creation 
and deletion, respectively. A schematic illustration of both pro¬ 
cesses is shown in Fig. in which we show how the periodic 
box moves along the orbit and the effect of the rotation scheme 
within the periodic box. 


2.1. Rotation scheme and the reference frame 


To simulate the tidal forces upon the dwarf galaxy, we apply the 
standard fictitious forces throughout the periodic box to simulate 
a solid rotation of the box (a non-solid rotation is also possible, 
although this introduces further issues, e.g. [Huber & Pfenniger[ 

[200T] ). 


To apply the fictitious forces, we apply the standard equa¬ 
tions of motion in a rotating reference frame, describing the ro¬ 
tation through the use of quaternion^ around a pivot point in¬ 
side the box. Although not commonly seen in astronomy, quater¬ 
nions are extensively used to describe rotations in programming. 
Quaternions, which are an extension of complex numbers, are 
able to describe 3D rotations with the useful properties that they 
can be easily updated with time and are reasonably forgiving 
of the gradual accumulation of errors. As the rotation involves 
renormalisation the errors only accumulate in the rotation di¬ 
rection/angle and do not change the magnitude of any vector. 
Throughout we write equations that involve the rotation of a vec¬ 
tor by a quaternion. Here when written in the same equation as a 
quaternion, a vector x should be understood to be the quaternion 
(x,0). 


2.1.1. Reference frames 

We locate the fixed pivot point within the periodic box at Xp = 
(Xp,yp,Zp), which represents the (moving) point relative to 
the host galaxy. For clarity, we use upper case variables when 
discussing things in the reference frame of the host galaxy and 
lower case within the periodic box. At the start of the simulation, 
this pivot point is set to be the position of the particle with the 

^ We write quaternions in the vector form q = (x, w) = (vi, ^ 2 , = 

xii + X 2 j + x^k + w, where f = f - ijk - w is the real 

part of the quaternion q and Xi,X2, V3 are real numbers corresponding 
to the magnitude of the imaginary components. The addition of two 
quaternions qi = (xi,wi) and q 2 = (X 2 ,W 2 ) is given by qi + q 2 = 
(xi + X 2 , Wi + W 2 ). The (non-commutative)multiplication of qi and q 2 
is given by qiq2 = {w\X2 + W2X1 + Xi X X2, W1W2 - Xi • X2). The inverse 
of a quaternion q = (x, w) is given by q"^ = (-x, K;)/||q|p, where ||q|| = 

yjx\ + x\ + x\ + uA is the norm of q. A rotation of a vector x by an 
angle 6 about an axis y (resulting in vector z), can be represented by 
(z, 0) = q(x, 0)q"\ where q = (y sin[^/2], cos[^/2]). 
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Fig. 1. A schematic illustration of the method. This is a representation of the simulation RPS-lOO described below ~400 Myr after infall. The left- 
hand panel shows the predicted orbit for 7 Gyr and the position of the rotating frame at the snapshot. The arrow shows the direction of the dwarf’s 
orbit. The right-hand panel is a quiver plot showing the position and velocity of a 0.5 kpc central slice of gas, with contours of the projected dark 
matter density shown in red. For conceptual ease, the velocity in an inertial frame is shown, which is related by a transformation to the rotating 
frame described in Eq. 


minimum potential inside the periodic box. As many equations 
arising from the rotation of the box depend upon the relative 
distance to the pivot, for clarity we use an asterisk to denote the 
difference between a vector and the pivot, e.g. X* = X - 
where necessary throughout the paper. The position of any other 
point, X, within the periodic box corresponds to the point X in 
the reference frame of the host galaxy by 

X* = q-'x*q, (1) 

where q is the relevant quaternion rotating a vector from the host 
galaxy frame to the periodic box. We choose q to be the unit 
quaternion that maps the unit velocity vector of the pivot to 
one of the axis in the periodic box, which we arbitrarily choose 
to be -y. 

2.1.2. Equations of motion 

Within the periodic box the velocity is calculated in the usual 
way according to the standard equations of motion for the rotat¬ 
ing reference frame 

V = V* = qV*q“^ - tu X X*, (2) 

where oj is the angular velocity within the rotating reference 
frame. Because of the dependence on the position, the veloc¬ 
ity of particles that cross any boundary need to be recalculated. 
This recalculation is not hugely important and, in practice, the 
gas pressure has a tendency to correct the particle velocity; gas 
that crosses a boundary is unlikely to interact with the dwarf, but 
correcting the gas in this way reduces the error associated with 
the periodic box. 

Each particle within the box experiences an acceleration 
force a, dependent upon the potential the potential of the host 


galaxy, O; the angular velocity; the acceleration, present on a 
particle inside the periodic box at a position x and velocity v, is 
given by 

a = a* = q[-VO(X)-A4q"‘-2wxv-wx(wxx*)-wxx*, (3) 

where O (X) is the potential at the point X, the acceleration of 
the pivot and the standard Coriolis, centrifugal, and Euler forces. 

The acceleration of the pivot, = -VO(Xp) -r A^ c. com¬ 
prises the gravitational acceleration due to the potential of the 
host galaxy plus an added correction term, A^^c, to account for 
the drag the dwarf experiences through the orbit. This drag and 
the effect on the dwarf’s orbit is discussed further in §4.6| 

The correction term is ad hoc, intended to disappear when 
not needed and be sufficiently small when needed to minimise 
numerical errors that occur when Ap^c is large (most of these er¬ 
rors manifest as issues with particle creation/deletion, discussed 
below). The correction manifests as a moving of the box off its 
original orbit following the new trajectory of the dwarf, modi¬ 
fied by ram pressure and tidal stripping. The correction term is 
calculated by finding the average position X 64 of the 64 stellar or 
dark matter particles that had the lowest potential at infall. The 
choice of 64 particles balances the risks between just taking the 
average of all particles, which results in a poor approximation 
when the galaxy becomes tidally disrupted, and choosing only a 
handful, which can focus on tidal tails if any of the particles are 
ejected from the core. In the simulations discussed here we use 

Ap,c = 0.54 Gyr-i + 0.022 GyrV (x*^, ¥* 4 ) X*^, (4) 

where f(x, v) is a function evaluated in each Cartesian dimension 
that returns 1 if, > 0 or |l>| <2 km s“\ and otherwise 0 . 
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2.1.3. Computing the rotation of the reference frame 

As q is not static, it must also be updated every time step. To do 
so, at every time step we first calculate the the angular velocity, 
n in the host galaxy frame, and oj = in the periodic box. 

The parameter O is calculated whenever the velocity is updated 
and is given by the standard relation 


Once the angular velocity is known, q is given by the useful 
relation 

q = O.SOq, (6) 

which allows the rotation to be updated each time step in time 
with the velocity. After each update in q it is renormalised, en¬ 
suring llqll = 1, and ensuring that errors only accumulate in the 
angle and do not artificially increase the magnitude of the vec¬ 
tors. 

That the rotation quaternion changes every time step and that 
the velocity of a particle appears in the calculation of the ac¬ 
celeration is problematic. Gadget-2 and by extension GEAR re¬ 
lies on a standard Kick-Drift-Kick integrator, which relies on the 
Hamiltonian of the system being separable into independent po¬ 
sition and momentum coordinates. Unfortunately, the standard 
Hamiltonian for this system breaks this condition for rotating 
reference frames and, as such, we lose the symplectic nature of 
the leapfrog integrator. In an effort to minimise the error result¬ 
ing from this, we make a slight alteration to the calculation of 
the acceleration when implemented in the code. 

Normally, the velocity is updated by an acceleration term, 
which only depends upon the position 

Vf+i/2 = ^t-1/2 + a(xOAr. (7) 

However, in the rotating frame it is clear from eqnj^that the ac¬ 
celeration term also includes a velocity dependence and terms 
dependent upon the changing rotation quaternion and angular 
velocity. We instead update the velocity through an iterated ap¬ 
proximation of the velocity at the time step, which performs rea¬ 
sonably well in simulations of the tidal force (tests of this method 
and others are provided in appendix]^. Care should be taken to 
note that this scheme is not the best performer at large step sizes, 
such as in dark matter only simulations. First, we assume that 
the change in rotation quaternion is negligible across any time 
step, secondly we assume that the acceleration in oj is slow com¬ 
pared to the time step size. These two assumptions hold well at 
time steps determined by the Courant condition inside the simu¬ 
lation. The first assumption leads to qr,pred = and the latter 
assumption means that we can approximate tu at a time t by 

t^f,pred = (8) 


2.2. Creation 

To simulate the increase in density over the orbit of the dwarf, 
particles are created and inserted into the box through the bottom 
side, we describe the exact mechanism of creation in appendix 

The velocity of the particle is then determined by the veloc¬ 
ity of the pivot point and the correction to the rotating frame of 
reference, 

y = (0, llVpIl, 0) - w X fx' - (0, i/p, 0)1. (12) 

As the orbit can change over the scale of the box this injection is 
not perfect though and the gas hits the dwarf at a slight angle. 

The remaining properties, temperature, metallicity, pressure, 
entropy, etc., are determined according to the properties of the 
chosen halo profile. Our choice of profile is a hot gaseous halo in 
hydrostatic equilibrium described in |3.2[ but in practice can be 
any function, with the caveat that a profile that rapid changes in 
pressure results in pressure waves. These errors and simple tests 
of the changing density profile are discussed further in Appendix 
1^ Notably, these errors depend on the sign of Rp and therefore 
partially cancel out, the results therefore may have an error in 
timing, but qualitatively are reduced over a full orbit because of 
this sign dependence. 

2.3. Deletion 

To prevent particles from reappearing in the periodic box and in¬ 
teracting with the dwarf, they are deleted when leaving the rear 
of the box, or in the case of dark matter and stellar particles, 
from the sides or front of the box. As gaseous particles possess 
a pressure, we consider them separately from dark matter and 
stellar particles when considering which particles to delete. Any 
gas particles (halo or dwarf), which upon moving drift within 
L/lOO of the rear box edge on that time step, are deleted. The 
choice of deleting them slightly before the edge of the box serves 
two purposes, one is to minimise the amount of neighbour par¬ 
ticles occurring across the rear boundary, and secondly, to pro¬ 
vide a small boost to particles about to leave the box. This slight 
vacuum does produce an increase in density waves of particles 
entering the box since they see the absence of particles on the 
other side, but these oscillations are greatly damped by the time 
the particles reach the dwarf. The lack of neighbours across the 
edges minimises errors that can occur near perigalacticon when 
the dwarf’s orbit changes rapidly (and so the particles may have 
minimal velocity in the y direction). 

Any dark matter or stellar particles are deleted if they drift 
within L/500 of any box edge or L/100 of the rear box edge. 
This extra deletion prevents a spurious build-up of particles at 
the edges of the box, which can occur when tidal tails pass the 
boundaries of the box. 

To assist in understanding, we illustrate this mechanism in 
Fig. in which we show which particles would be deleted in 
future time steps if they were gas, dark matter, or stellar particles. 


where tUoid is the change in oj at the previous time step. The 
velocity is then updated from a time t-l 12 At to a time t+\ 12 At, 
according to 

Vf,pred = v^_i /2 + a(q^_i/ 2 , X^, V^_i/ 2 , tu^_i/2)Ar/2, (9) 

Vr,pred ~ V^-1/2 + a(q^_l/2, X^, V^^pj-e(5, tU^^pj-ed)^^/2, (10) 

^t+l/2 — V^-1/2 + a(q^_l/2, X^, V^^pred? tU^,pred)^^- (H) 
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Fig. 2. Schematic of particle deletion. Shown in a 50 x 50 kpc box. 
The deletion regions (exaggerated by a factor of 10 from that used) are 
shown with a red dashed line (for deletion of all particles) and green 
dotted lines (deletion of stellar/dark matter particles). Particle positions 
and velocities are shown with arrows. Particles that would be deleted 
on the next time step if they were stellar or dark matter particles are 
shown in green. Particles that would be deleted on the next time step 
regardless of their type (gas, stellar, or dark matter) are shown in red. 
The remaining particles, which are not deleted regardless of type, are 
shown in black. 


3. Simulations 


To investigate the impact of ram pressure and tidal forces upon a 
dwarf galaxy, we investigate two Sextans-like dwarf galaxies, 
one with a pseudo-isothermal dark matter profile reminiscent 
of today’s dwarf ( [Battaglia et al.|[20TT] ) and one with a denser 
NFW profile. In both cases, we use the same pseudo-isothermal 
gas profile. We investigate three scenarios, the dwarf in isolation 
and as it undergoes ram pressure throughout two different orbits, 
both of which have an apogalacticon of 200 kpc and differing 
perigalacticons of 50 kpc and 100 kpc. 

The dark matter and gas profile parameters (see Table [T]) and 
parameters for star formation and feedback are calibrated such 
that this isolation model produces a dwarf with chemodynami- 
cal properties similar to Sextans dSph after 6 Gyr total. Based on 
this calibration throughout all models, we use a star formation ef¬ 
ficiency of c* = 0.01, a supernova energy of 6sn = 0.5 x 10^^ erg, 
a density threshold for star formation of pth = 0-01 cm“^ , and 
temperature threshold of Tth = 3 x 10"^ K. Changing the den¬ 
sity or temperature thresholds has minimal impact on the fi¬ 
nal results at the resolutions considered here (see e.g. [Revaz & 


[Jablonka[|MT2| . As per [Revaz & Jabio5S[ ( |20T^ and [NRJ14| 


below T = 10^ K, we use the metallicity-dependent cooling ac¬ 
cording to [Sutherland & Dopita[ ( [1993] ). Below this, we also ac¬ 
count for H 2 and HD cooling and cooling due to fine structure 
line transitions of oxygen, carbon, silicon, and iron, according to 
IMaio et al.| ( |2007l l. 

The initial mass function (IMF) is chosen to match that of 
[Kroupa[ ( [2001[ ), and we carefully track the number of Type la 


and Type II supernovae that explode at each time step calculat- 
ing the rates according to the random scheme of [Revaz et al. 
( 2015[ ). This treatment has shown good energy conservation and 
works well over a large variety of masses reproducin g accurate 
stellar populations when compared to observations ( [Revaz & 
Jablonka|2012[ [Revaz et al.|2015[ ). Metals are diffused through¬ 
out the nea rby gas with a smooth metallicity scheme ( [Wiersma] 
et al.|2009) ) whereby metal contents are calculated via the SPH 
scheme at any individual time step, which reduces the dispersion 
in [of/Fe] at low metallicity. We show the resulting star formation 
history and metallicity profiles in Fig.|^ 


3.1. Initial density profiles 

In half of our models, the dwarf’s dark matter profile is initially 
described by a pseudo-isothermal profile, 


p(r) = 


Pc 


1 + r'^jrl ’ 


(13) 


where pc is the core density and is the core radi us. This choice 
of profile has much observational support (e.g. [de Blok et al.[ 
[2008[ [Walker & Penarrubia|2011 ) although if the true profile of 
dwarf spheroidals is clos er to NFW (which can not be ruled out; 
[Breddels & Helmi|[2013[), it is likely to overestimate the effects 
of tidal interactions ( [Kazantzidis et ar][2013[ ). To examine the 
effects of the profile choice, we also consider an NFW profile 
( [Navarro et al.|1997] ) for the dwarf’s dark matter, 


p{r) = ps 


r/vs 


(1 -r rlvsfrlrs ’ 


(14) 


where ps is a scale density and a scale radius. In both cases, the 
velocity of the dark matter is determined from the Jeans equa¬ 
tions. 

The gas associated with the dwarf is initially set at the same 
position as the dark matter and slightly rotated to avoid overlap¬ 
ping particles, with a temperature determined by converting all 
the kinetic energy. The kinetic energy is determined by the Jeans 
equations into thermal energy, thereby giving the gas zero veloc¬ 
ity initially. We show the simulation parameters in Table and 
the model differences in Table |2| 

To test convergence, we also ran all simulations at higher 
resolutions with particle masses a factor of four lower. We find 
very little difference dynamically, with gas being removed at the 
same time and the final state of the dwarfs being dynamically 
similar. However, as increased resolution can cause a spurious 
increase in the scatter of [tr/Fe] ( [Revaz et al.|2015[ ) because of 
scatter in the IMF, we use the larger particle gas mass. 


3.2. Host galaxy 

The host galaxy is described by a multi-component static po¬ 
tential within which a dwarf experiences tidal effects and a hot 
gas halo, which provides ram pressure stripping. The potential of 
the host galaxy consists of a bulge and disk modelled as Plum¬ 
mer profiles with parameters chosen to match the Milky Way 
rotation curve of [Xue et ni(|2008[ ) and an NFW dark matter halo 
with a virial mass of 8x10" M© and concentration c = 21 as per 
Kafle et al.[ ( |2QT^ . We model the disk as a spherical potential, 
because large distances are involved (at least 50 kpc) and this 
approximation makes little difference to the final form. The po¬ 
tentials of the bulge and disk have minimal effect for most of the 
orbit and could be left off, however, the inclusion of their mass 
does result in a slight increase in velocity at perigalacticon and 
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Pseudo-isothermal 



NFW 



Fig. 3. Properties of the pseudo-isothermal dwarf (top) and NFW dwarf (bottom) in isolation after 6 Gyr of evolution. The properties of the Sextans 
dSph are shown with a red dashed line in the [Fe/H] histogram (mid dle panel) ([Battaglia et al.|2011|). The [Mg/Fe] v [Fe/H] of available RGB 
stars ( [Shetrone et al.|200T]|Aoki et al.|2009[|Tafelmeyer et al.|2010 JKirby et al.poTT |Theler et al.|20T5] ), selected as per |Theler et al.| ( [20T5] ), 


are shown with points in the third panel. The Mg values are derived from one line and show much larger scatter at low a that is not observed (see 
|Theler et al.|2015) . The median [Mg/Fe] in each [Fe/H] bin of the stellar particles is shown with a cyan line in the same panel. This choice of 
parameters achieves the approximate metallicity distribution (although missing the most metal rich stars) and the turnover in [cr/Fe] of Sextans 
dSph. 


Table 3. Host halo parameters 



(p(R) Form 

Parameters 

Reference 

Bulge 

-GM/ V7?2 + a2 

M= 1.3x 10^0 Mo 
a = 0.5 kpc 

[1] 

Disk 

-GM/ + ^2 

M = 5.8 X 10^0 Mo 
a = 5 kpc 

Mvk = 8 X 10" Mo 

[1] 

Halo 

-GMvir ln(l + 7?c//?vir)/[ln(l + c) - c/(l + c)] 

c = 21 

Ryir = 240 kpc 

[2] 


Notes. References: [l] |Xue et al.| ( [20081 ) [2][Kafle;e^([2^ 


because of the strong dependence of ram pressure on the veloc¬ 
ity, we include them for the sake of completeness. The forms and 
parameters of which are given in Table 


The hot halo of the galaxy is considered to be an isothermal 
halo of 2 X 10^ K in hydrostatic equilibrium with the potential. 
This choice only requires a scale density, which we set by de¬ 
termining that the electron number density at 50 kpc is given by 
He = 2x 10~^ cm“^. The profile is shown in Fig. and compared 
with the profile of Miller & Bregman] ( |2013[ red dashed line. 


with p errors in red dotted line) and |Gatto et al. (|2013|). Th is pro¬ 
file gives a high density compared to Miller & Bregman] ( [201 3| ) 
profile, particularly at high galactocentric radius ( R^lOO kpc). 
This profile is approximately the same as that of jGatto et al. 


( |2013| in which the blue points) are necessary to strip the Sex¬ 
tans and Car ina in one perigalactic passage , and is slightly below 
that used by [Bland-Hawthom et al.j pOOVj ) to reproduce the Ha 
emission along the Magellanic Stream. 
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Table 1. Initial density parameters 


Parameter 

Value 

4^tot 

3.5 X 10« Mg 

4^DM,tot 

2.98 X 10^ Mo 

4^gas,tot 

0.53 X 10^ Mo 

4^DM,part 

39.5 X 10^ Mo 

^gas,part 

6.0 X 10^ Mo 

r c,gas 

2 kpc 

Pc,gdLS 

5.4 X 10^ Mo kpc-3 

Pseudo-isothermal models 

rc,DM 

2 kpc 

Pc,DM 

3.0 X 10® Mo kpc-3 

NFW models 

G,dm 

1 kpc 

P5,dm 

2.3 X 10’ Mo kpc-3 


Notes. Simulation parameters for the dwarf galaxy. The parameter Mtot 
is the total mass of the system (excluding the low-density medium), 
^RM,tot the total dark matter mass, Mgas,tot is the total gas mass (exclud¬ 
ing the low-density medium), A^part is the total particle number (exclud¬ 
ing the low-density medium), MDM,part is the dark matter particle mass, 
^gas,part is the gas particle mass (of all gaseous particles), is the star 
formation parameter, 6 sn is the energy released per supernova, r^^gas is 
the core radius of the gas in the dwarf, pc,gas is the core gas density of 
the dwarf. 


Table 2. Simulation parameters 


Model 

Profile type 

Perigalacticon 

^peri 

[kpc] 

Apogalacticon 

^apo 

[kpc] 

ISO 

Pseudo 

- 

- 

RPS-50 

Pseudo 

200 

50 

RPS-100 

Pseudo 

200 

100 

NISO 

NFW 

- 

- 

NRPS-50 

NFW 

200 

50 

NRPS-100 

NFW 

200 

100 


Notes. Model peri- and apogalacticons and whether they have a pseudo- 
isothermal or NFW dark matter profile. 


3.3. Model Initialisation 

For the first two gigayears, all simulations are identical and 
evolved without any external potential but with a surrounding 
warm medium. The dwarf is evolved in the centre of a periodic 
box of side length L = 50 kpc. Outside the dwarf a warm (10"^ K) 
low-density (3.7 x 10“^ cm“^) medium fills the box, this den¬ 
sity is equal to the density of the hot halo at 200 kpc but of 
a much cooler temperature and therefore lower pressure. This 
kind of medium is used to simulate the effect of the intergalactic 
medium. The isolation model is then evolved for a further 12 Gyr 
in this state. 

3.4. Insertion into the hot halo 

For the RPS models, after 2 Gyr the dwarfs are extracted from 
this periodic box by selecting only particles within a radius 
of 10 kpc from the centre. Only a handful of particles origi¬ 
nally associated with the dwarf are beyond this radius, and the 
exact choice does not impact the results appreciably. This re¬ 
moval of particles outside this radius minimises the slight over- 
pressurisation of the hot halo that would result if the particles 



Fig. 4. Halo density profile for our chosen model (black line). The best- 
fit p profile of |Miller & Bregm^ ( |2013| is shown with a red dashed line, 
with dotted lines indicating the error in p. The density of the halo deter¬ 
mined by the simulations of |Gatto et al.| ( |2013] ) for Carina and Sextans 
to be stripped in one passage are shown with blue points. 

are kept. The dwarf is then reinserted in the centre of a new pe¬ 
riodic box of the same size L = 50 kpc filled with a hot gaseous 
medium with a relative velocity according to the halo density 
profile at R = 200 kpc (n^ = 3.7 x lO""^ cm“^) and correspond¬ 
ing orbital velocity. We do not remove any of this medium when 
inserting the dwarf. 

Not removing this medium produces an overlap between 
the hot medium and the dwarf, however, at the slow (but non- 
negligible) velocity the dwarf has at apogalacticon, and because 
of the hot medium having extremely inefficient cooling, the 
dwarf quickly reaches pressure equilibrium with the external 
medium without losing gas or accreting the hot medium. 

4. Results 

As the dwarf undergoes extended evolution inside a hot halo, its 
chemodynamical properties post-infah are dramatically altered. 
The gaseous profile is affected almost immediately, the dark mat¬ 
ter profile is altered as gas is evacuated, and tides take place. 
These profile changes also alter the star formation and chemical 
evolution of the dwarf, altering the properties that would be ob¬ 
served today. Drag as a result of the ram pressure also changes 
the orbit of the result, reducing the perigalacticon and apogalac¬ 
ticon of the dwarfs, suggesting that dwarfs may arrive slightly 

later than pure kinematic studies may suggest (e.g. IRocha et al.l 

- 

4.1. Centre calculation 

For every property except the global star formation (shown in 
|4.2| ), we calculate the centre of the dwarf by applying a Gaus¬ 
sian smoothing function to the three-dimensional density his¬ 
togram of 64 pixels (originally over the box) with a standard 
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deviation of the Gaussian of 6 pixels. We then refine around the 
densest point, halving the size of each dimension considered in 
each iteration until a region of (6 kpc)^ is considered, and then 
perform the same operation again. From this densest point we 
do a simple mass weighted average position over all particles 
within 3 kpc of this point. This procedure has been found to be 
relatively stable around disrupting dwarfs in |NRJ14] and suc¬ 
cessfully avoids dense points in tidal tails in favour of the centre 
of the dwarf. 


4.2. Star formation 

Upon entering the confining medium of the halo where tem¬ 
perature greatly exceeds that used to simulate the intergalactic 
medium, star formation increases dramatically as a consequence 
of the increase in pressure. It has been suggested that this in¬ 
crease in star formation happens in the Milky Way satellites, 
based on comparisons between kinematic properties and star for¬ 
mation histories ( [Rocha et al.||2012| ). The total stellar mass is 
shown in Fig. which includes all stellar particles within the 
box at any time, including any that are formed in gas clouds 
as they are stripped from the dwarf. Here any drop in the stel¬ 
lar mass represents the tidal stripping of stars from the dwarf, 
which can take place even while star formation continues within 
the box. 

For the pseudo-isothermal models, the isolation case, ISO, 
continues a steady approximate continuous increase in stellar 
mass. Post-infall, the RPS-50 model begins to experience a large 
jump in star formation mass, before rapidly dropping just after 
the first perigalacticon passage (a sign of extensive tidal strip¬ 
ping as stellar particles are removed from the simulation). This 
bursty behaviour is refiected in the star formation rate shown in 
Fig.g 

The RPS-100 model continues star formation for a much 
longer period, but undergoes two separate phases. Firstly, the 
star formation is enhanced similar to that of RPS-50, but after 
the first perigalacticon passage, a slight decrease is seen corre¬ 
sponding to the change in gradient in Fig. before the dwarf is 
tidally disrupted after the second perigalacticon passage. 

For the NFW models, because of the denser dark matter 
profile, the stellar mass tends to be higher than in the pseudo- 
isothermal models. Again, the isolation case, NISO, continues 
a steady approximately continuous increase in stellar mass. For 
those post-infall, however, the star formation is much more con¬ 
tinuous than in the pseudo-isothermal models. Still enhanced, 
the models continue star formation until in the case of NRPS-50 
it is suddenly halted, and the rate of star formation in NRPS-100 
drops after each perigalacticon passage. 

Here, the star formation is enhanced in NRPS-50 before be¬ 
ing shut down around 6 Gyr as gas is finally removed from 
the dwarf towards near the second perigalacticon (but see also 
|4.3.2| ). In NRPS-100, the enhancement is again seen, but after 
perigalacticon the star formation rate decreases slightly, corre¬ 
sponding to the change in gradient in Fig.|^ Before finally being 
stopped only a few hundred Myr before the end of the simula¬ 
tion. That a rapid change in star formation at perigalacticon may 
allow the determination of the time of perigalacticon passage, 
particularly for recently accreted dwarfs where the time resolu¬ 
tion of the star formation history is higher. 
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4.3. Mass evolution over time 

At the same time that the gas is moving within the dwarf, the 
mass contained within the central regions of the dwarf changes 
dramatically. Figs[7]shows the mass inside the central 1 and 3 kpc 
for the pseudo-isothermal and NFW models versus time. 


4.3.1. Dark matter mass 


The dark matter mass is stable (after initial settling) in the iso¬ 
lation model shown in panel (a), with slight oscillations in the 
mass as gas is blown out and reaccretes. As tidal forces begin to 
affect the dwarf in RPS-50 and RPS-100 near the first perigalac¬ 
ticon, a dramatic drop in dark matter mass is seen. The drop in 
mass at the same tim e as baryons is lost is not surprising hav¬ 
ing been explored in [Arraki et al.j ( |2014| ), who found that the 
removal of baryons from a dwarf tends to drop the density (and 
therefore mass) by a factor up to (1 -r /b)"^. The effects here how¬ 
ever are much stronger than the -50% mass loss predicted by 
Arraki et^ for total baryon loss, and may be the consequence 
of all the mass being lost near perigalacticon. 

This reduction is enormous when considering the closest 
either model gets is 50 kpc, a distance far from sufficient to 
tidally destroy dwarfs when ram pressure stripping is not in¬ 
cluded ( |NRJ14| ). Here the RPS-50 model is effectively destroyed 
after one passage, and the RPS-100 model is reduced by almost 
an order of magnitude and exists as only a stream, as indicated 
by the dark matter mass being near zero within the central 3 kpc. 

The synergy between baryon loss and dark matter loss has 
also been noted in [Zolotov et al.j ( [2012[ ), where dwarfs with 
baryons lost a greater mass fraction than those in a simulation 
that only includes dark matter, particularly in those cases where 
dwarfs have a baryonic disk. Here, with no disk present we see 
similar synergy, and this drop in density occurs due to the alter¬ 
ation in the potential, resulting from gas being lifted to higher 
orbits from superno va feedback ( [Govemato et aU2()T0[ [Pontzen 
& Govemato 2012) and ram pres sure stripping That this tidal 
effect is so large, completely removes the pseudo-isothermal 
dwarfs resemblance to a classical dwarf spheroidal today, sug¬ 
gesting that Sextans at least looked remarkably different at infall 
to its present day structure. 

Compared to ISO, NISO, shown in panel (b), shows a higher 
mass of dark matter in the central 3 kpc, but the same fluctua¬ 
tions, corresponding to bursts of star formation and outflows of 
gas, are seen as in the ISO case. As with the pseudo-isothermal 
models, dramatic drops in dark matter mass are seen when the 
dwarfs approach perigalacticon, in addition to the normal 1 kpc 
fluctuations. Because of the denser inner region of the dwarf, 
these dwarfs are more resistant to tidal stripping and, in both 
cases, take almost twice as long to be transformed into streams. 

To ensure the mass loss is a consequence of the ram pressure 
stripping and not a numerical artefact we have also simulated the 
evolution of a baryon free dark matter dwarf and a dark matter 
dwarf moving through a baryon free halo and therefore not expe¬ 
riencing ram pressure. In both cases, these dwarfs survive with 
a noticeable dark matter core until the present day, suggesting 


^ Although the mechanisms from supernova feedback are likely to 
lower the central density independently, in these simulations they do 
not occur rapidly enough to result in a net change. However even with 
an adiabatic change, a core can be created because the gas is prevented 
from reaccreting by ram pressure stripping that prevents the dark mat¬ 
ter from adiabatically contracting and consequently the density is still 
reduced. 
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(a) Pseudo-isothermal (b) NFW 


Fig. 5. Mass of stellar particles inside the simulations at any time (cumulative star formation minus stellar mass lost to tides and supernova 
explosions). ISO is the pseudo-isothermal model in isolation, RPS-50 with a perigalacticon of 50 kpc, and RPS-lOO with a perigalacticon of 
100 kpc. NISO is the NFW model in isolation, NRPS-50 with a perigalacticon of 50 kpc, and NRPS-100 with a perigalacticon of 100 kpc. 
The time of infall (before which all simulations are identical) is shown with a black dashed line, approximate perigalacticons are shown with a 
dot-dash line in the relevant colour and apogalacticons in a dotted line. In all cases, the peri-/apogalacticon shown is of the pivot and the actual 
peri-/apogalacticons of the dwarf may differ slightly in time. 


that ram pressure exacerbates tidal stripping. We describe these 
simulations in Appendix [D| 


4.3.2. Gas mass 

In contrast to the relatively stable nature of the dark matter in 
ISO, the gas mass in all pseudo-isothermal simulations, shown 
in panel (c), changes dramatically with bursts of star formation, 
particularly in the central 1 kpc, which changes with star forma¬ 
tion by a factor of two. In the ram pressure simulations, the cen¬ 
tral 3 kpc also experiences dramatic changes of 30-40%, with 
the central 1 kpc experiencing changes by more than a factor of 
three. Although only a few bursts take place, a general trend is 
seen, with slightly less material reaccreted each time, particu¬ 
larly to the inner 1 kpc before finally being stripped near the first 
perigalacticon. 

In the RPS-50 model, a spike is seen around 5 Gyr (3 Gyr 
post infall), which is caused by stripped gas re-entering via the 
side of the box and within the central 3 kpc of the dwarf. As this 
gas has a large velocity, it never forms stars that remain bound to 
the dwarf, and subsequently does not impact any future analysis. 
Around 8 Gyr (6 Gyr post infall) and 11 Gyr (9 Gyr post infall) 
a slight increase in the mass is seen because the dense galactic 
corona has enough mass within the central 3 kpc to be noticeable 
and the rise in gas mass at this time is not an accretion event. 

For the NFW models, the gas mass, shown in panel (d), 
shows similar trends as the pseudo isothermal models, Once 
again, the compressive force increases the star formation and the 
oscillations of this gas. In the NRPS-50 model, a clear drop in 
the gas mass is seen at perigalacticon, soon after a large amount 


of gas is then reaccreted at apogalacticon (~6 Gyr). This gas is 
held slightly off centre throughout the dwarf’s passage before the 
lower ram pressure forces allow it to recall into the central po¬ 
tential well of the dwarf, a mechanism for star formation bursts 
suggested previously in [Nichols et al.| ( ~2012| ). This is the only 
model (including those not shown in this paper, with variations 
in star formation efficiency, supernova feedback, initial condi¬ 
tions etc.) in which this kind of mechanism appeared, and there¬ 
fore any reaccretion event like this is likely to be rare. The small 
increases seen at the perigalacticon passages are also visible. The 
NRPS-100 model has a comparatively simple history, gradually 
decreasing in mass as ram pressure strips the outer edges of the 
gas until the dwarf is tidally destroyed and all gas is loss. The 
complicated nature of the gas removal suggests that any simple 
method to remove gas in semi-analytic simulations will likely 
end up as a poor approximation. 


4.3.3. Stellar mass 

The stellar mass in the pseudo-isothermal models, shown in 
pane l (e), increases with bursts of star formation, discussed in 
|4.2| with some change within the central 1 kpc as gas expul¬ 
sion drags some of the stellar mass to higher orbits, before the 
reaccretion brings material back. The stellar mass in the central 
3 kpc for the ISO model grows monotonically with time, with the 
gas expulsion failing to drag any stars further out than this. As 
the ram pressure models approach perigalacticon, a large stellar 
mass is lost from the inner central kiloparsec. For RPS-50, these 
stars become unbound as the dwarf is rapidly transformed into a 
stream, but for the other ram pressure models most of these stars 
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(a) Pseudo-isothermal 

Fig. 6. The star formation rate inside the simulations. Lines and colours are the same as in Fig.[^ 


(b) NFW 


move to higher orbits outside the inner kiloparsec, with the stel¬ 
lar mass inside only suffering a small reduction until the dwarf 
is completely destroyed. 

Because of the denser dark matter profile, the stellar mass 
within 1 and 3 kpc in the NFW dwarfs, shown in panel (f), 
tends to be higher than in pseudo-isothermal dwarfs. As tidal 
forces begin to act on the dwarf, many stars move to higher or¬ 
bits within the dwarf, reducing the mass within 1 kpc even while 
the total stellar mass within 3 kpc continues to increase. This is 
partially counteracted in NRPS-50 by the reaccretion event drag¬ 
ging some stars back to lower orbits, in addition to the formation 
of new stars. As in the pseudo-isothermal models, the dwarfs 
eventually become streams and the stellar mass within the cen¬ 
tral 3 kpc drops to ~0 as a consequence of this. 

Ultimately all pseudo-isothermal and NFW models become 
tidal streams, however, the NFW profiles survive for a much 
longer time, surviving past the first perigalacticon with a large 
amount of gas. 


4.4. Spatial gaseous and stellar distribution 

As gas is stripped from the dwarf, it goes through numerous 
forms, which is shown for the pseudo-isothermal dwarf in Fig.[^ 
and for the NFW dwarfs in Fig. In Fig. 9, we show gas with 
a temperature below 10^ K through a 0.5 kpc slice of the dwarf. 
As would be expected, the loosely bound outer gas of the dwarf 
is quickly removed within a few hundred megayears. The in¬ 
crease in star formation and related feedback, resulting from the 
confining pressure, begins to blow out gas quite quickly. How¬ 
ever, because of the relative velocity of the medium, this gas is 
not reaccreted in spherical shells, and instead forms a fiattened 
bullet-like shape (see e.g. RPS-lOO at 3.5 Gyr or NRPS-100 at 
2.75 Gyr), with the leading edge of gas compressed by the am¬ 
bient medium. By 1.5 Gyr after infall (3.5 Gyr of evolution), the 
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dark matter is warped by tides and the removal of gas, becoming 
elongated in the direction of motion. At the same time, the gas 
is no longer centred on the dwarf, but is offset to the side that is 
opposite the direction of motion. 

This off-centre behaviour results in a noticeable change in 
the location of star formation (discussed below) and at the same 
time slows down the dwarf altering its orbit (this interesting fea¬ 
ture is discussed later in ^4.6| ). 


In the NFW models, a similar evolution is seen, except that 
the central gas is denser, and fails to experience the massive 
blowouts. Similar to the pseudo-isothermal models, slight dis¬ 
tortions in the outer dark matter profile are seen, and the gas is 
held slightly off-centre. As the dwarf begins to approach peri¬ 
galacticon, all the gas in RPS-50 is removed (last row), while in 
NRPS-50, RPS-lOO, and NRPS-100, some gas is able to survive 
and is held off-centre by the dwarf. 

The fact that the gas is no longer distributed in a roughly 
spherical fashion has a large impact on the positions of star 
formation. In Fig. we show the luminosity for the pseudo- 


isothermal models at the same times as the gas density in Fig. 

As the luminosity is dominated by young stellar populations, 
this traces the star-forming regions within the dwarf. As always, 
star formation takes place in regions where dense gas resides, in 
dwarfs undergoing ram pressure, which has the effect of creating 
obvious distortions from sphericity in the luminosity profile, par¬ 
ticularly in the pseudo-isothermal models. However, the dwarf’s 
gravity and tidal effects tend to hide these distortions with time. 


In the NFW models, shown in Fig. [TT] this effect is much 
more reduced with nearly all the stars held in the central region, 
with only a low brightness halo of star formation, as would be 
expected from the reduced blowouts. These blowouts still result 
in a slight distortion in the NRPS-100 model as in RPS-lOO, and 
similarly, the tidal tails of RPS-50 also appear inside NRPS-50. 
Of particular note is that as the last gas is expelled, stars continue 
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Fig. 7. Mass of dark matter (top row), stars (middle row), and gas (bottom row) within the central 1 and 3 kpc of the pseudo-isothermal and NFW 
dwarfs. The ISO model is shown with a black line, RPS-50 with a blue line, and RPS-lOO with a red line. The NISO model is shown with a green 
line, NRPS-50 with a cyan line, and NRPS-100 with a magenta line. The 3 kpc line is shown with thick lines and 1 kpc with thin lines. The infall 
time is shown with a vertical black dashed line. 
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"2.00 Gyr ISO RPS-50 RPS-lOO" 




' 2.75 Gyr 



Density (cm“^) 

Fig. 8. Density of gas (image) and dark matter (contours) at 750 Myr intervals for pseudo-isothermal models from infall; beginning at infall in a 
0.5 kpc slice through the centre of the dwarf. Only gas below 10^ K is shown. Each box is 15 kpc per side. The ISO model is shown in the left 
column, RPS-50 in the middle column, and RPS-100 in the right column. The gas density is listed on the colour bar and the dark matter contours 
(from thick to light) are 10^-^, 10^, and 10^-^ M© kpc"^. 
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Fig. 9. Density of gas (image) and dark matter (contours) at 750 Myr intervals for NFW models from infall until around the first perigalacticon, 
beginning at infall in a 0.5 kpc slice through the centre of the dwarf. Only gas below 10^ K is shown. Each box is 15 kpc per side. The NISO model 
is shown in the left column, NRPS-50 in the middle column and NRPS-100 in the right column. The gas density is listed on the colour bar and the 
dark matter contours (from thick to light) are 10^^, 10^, and 10^ ^ M© kpc“^. 
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Fig. 10. Luminosity of stars at 7 50 Myr intervals fro m infall for the pseudo-isothermal models. The stellar particles are smoothed to encapsulate 
the 16 nearest neighbours, as per Torrey et al. < |2015| ). Cool gas (T <10^ K) column density contours are shown in red. 


to form in the clouds, resulting in a (temporary) bridge of stellar 4.5. Metallicity 
material outside the dwarf, particularly visible in the last row of 

Fig. [TT] for the NRPS-50 model. The increase in star formation that takes place upon infall would 

normally suggest a quick evolution of metallicity as generations 
of stars pollute the interstellar medium. However, when a dwarf 
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Fig. 11. Luminos ity of stars at 750 M yr intervals from infall for the NFW. The stellar particles are smoothed to encapsulate the 16 nearest 
neighbours, as per Torrey et al. ( |2015| . Cool gas {T < 10^ K) column density contours are shown in red. 


falls into a hot halo, its metallicity no longer necessarily under¬ 
goes this rapid evolution due to star formation taking place in 
low-metallicity gas originally from the exteriors of the dwarf 
(discussed further in |4.5.1| ). We show the [Fe/H] in Fig. [T^ 


with the pseudo-isothermal models in the left panel and NFW 
models on the right. For each model, we consider the metallicity 
at the point of maximum stellar mass. 
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While the ISO and RPS-100 models both peak around 
IFe/H]—1.7, the RPS-50 model peaks at IFe/H]—2 despite 
both the RPS-50 and RPS-100 models possessing a greater stel¬ 
lar mass than the isolation model and star formation taking place 
over a number of gigayears. 

For the NFW isolation model, NISO, the IFe/H] peaks at 
higher metallicity IFe/H] —1.3, a consequence of the greater 
star formation. However once again the ram pressure model with 
a perigalacticon of 50 kpc peaks at IFe/H]—2 and at the time 
of the reaccretion event (slightly later than the peak of stellar 
mass but representing a more evolved dwarf) has a peak of only 
IFe/H]—1.6, which is exaggerated by the loss of stars due to 
tidal stripping. With its much higher star formation, the NRPS- 
100 model is able to evolve further in metallicity, reaching a peak 
of IFe/H] —1.1, and forms a significant fraction of stars above 
tFe/H]= -1. 

This halting in IFe/H] evolution is also seen in the IMg/Fe] 
evolution shown in Fig. Because of the increase in star for¬ 
mation, Type II supernova continue to enrich the medium in al¬ 
pha elements, resulting in a higher IMg/Fe] past the knee. This 
feature is common between the pseudo-isothermal and NFW 
ram pressure models and results in the isolation models having a 
larger fraction of low alpha stars (IMg/Fe] < 0) than any of the 
ram pressure models regardless of the IFe/H] distribution. 


4.5.1. Star-forming gas 

Although only a minor component (particularly for the NFW 
models), it can be seen in Fig. [T^ that stars of low-metallicity 
IFe/H] < -3 continue to form after infall. As all the central gas 
on infall possess a metallicity above this, at least some stars must 
be formed from gas originating in the outskirts of the pseudo- 
isothermal dwarf, and this is a likely cause for the halting in 
metallicity evolution. In Fig. we show the radius of gas at 
2 Gyr that goes on to form stars in the future. It is immediately 
clear that nearly all gas that goes on to form stars in the isolation 
case is located within the central 2.5 kpc of the dwarf. How¬ 
ever, for the RPS-50 and RPS-lOO model -25% of the stellar 
population is formed from gas outside this region. For the NFW 
models, the formation of stars of low metallicity is less com¬ 
mon, with only a few low-metallicity stars formed post infall. 
Once again, the radius of gas that goes on to form stars is higher 
for the ram pressure models than the isolation model. Here, the 
NRPS-50 model is similar to the NRPS-100 model despite the 
reaccretion event that occurs, which is indicative of most of the 
gas at large radii already being lost before perigalacticon. 

This gas from the outskirts appears to form gas through two 
different mechanisms. The primary mechanism is gas in the di¬ 
rection of travel is compressed, increasing in density and allow¬ 
ing it to eventually form stars. A secondary mechanism, how¬ 
ever, is that because of the dwarf’s gravity field, gas that is ex¬ 
pelled towards the front is perturbed enough to collide at the rear 
of the dwarf, allowing the gas to be reaccreted at the rear, effec¬ 
tively forming Bondi-Hoyle accretion. 


4.6. Orbital changes 

As evident from |4.3.2[ the stripping of the gas from the dwarf 
is not instantaneous, and for a large period of time, the gas is 
held off-centre as the dwarf’s gravity balances the forces result¬ 
ing from the ram pressure (see in particular the third row of Fig.[^ 
for an example of this). This gravitational restoring force results 
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in the gas of the dwarf acting like a parachute slowing down the 
dark matter and stellar content of the dwarf over many gigayears. 

The net impact of this off-centre gas is in the rotating frame 
of reference, a transfer of momentum and energy from the hot 
medium to the gas (gas to hot medium in the host galaxy refer¬ 
ence frame), and then to the dark matter and stellar component 
of the dwarf through the gravitational restoring force, allowing 
the gas of the dwarf to act as a parachute slowing the entire dwarf 
down. This transfer of energy results in a compression and heat¬ 
ing of the hot medium (see also [Roediger et al.|[MT5] ), but of 
larger impa ct here, a change in the orbit of the dwarf an effect 
also seen in [Smith et al.|p01^ . Because of the large mass con¬ 
trast between the dwarf and the galaxy and the static nature of 
the potential used, this effect is purely due to the ram pressure 
force. We verify that this effect is not due to the rotation scheme 
or the creation/deletion scheme in Appendix by performing 
two tests: firstly a dwarf in the hot coronae of the host halo, but 
without any baryons of its own; and secondly a dwarf in the ro¬ 
tation scheme without this hot coronae, which experiences only 
tidal forces. In both cases, the dwarf’s orbit remains constant 
with time, suggesting the effect here is due to the interaction be¬ 
tween the hot coronae and the gas of the dwarf. 

We show the change in specific energy (kinetic plus potential 
energy) of the centre of the dwarf and how it affects the orbit of 
a dwarf as a function of time in Fig.[^ The energy is calculated 
by examining the velocity and potential energy of the centre of 
the dwarf, excluding rapid changes due to the jumps in the centre 
determination, and is smoothed with a window function over 20 
snapshots (94 Myr). Selecting the centre of the dwarf means we 
do not have to worry about particle deletion over time, but this 
choice increases the noise at which the part of the dwarf that is 
considered the centre changes. Here, a general downwards trend 
is visible, as the dwarf loses energy and becomes more bound to 
the host halo. 

For RPS-50 and NRPS-50, which lose all their gas at the 
first perigalacticon passage, this loss of energy mostly occurs 
just prior to this passage, suggesting the dominant effect is re¬ 
ducing the apogalacticon. The longer period over which gas is 
removed from RPS-lOO and NRPS-100 is seen as the continual 
drop in energy over ~8 Gyr. As the dwarfs become streams, the 
calculation of the centre becomes quite uncertain. Resulting in 
spurious jumps in energy, these final jumps greatly exaggerate 
the (slight) gain in energy that can occur in the centre due to 
reaccretion of gas after a burst of star formation. A slight gain is 
seen as reaccreted gas was preferentially pushed in the direction 
of motion. 

At each snapshot, we also calculate the peri- and apogalac¬ 
ticon of the orbit of the dwarf, based on its current position and 
energy. Here the perigalacticon is slightly decreased, however, 
the apogalacticon experiences a -10%—25% drop. This is par¬ 
ticularly evident in the NRPS-100 model, where the extended 
loss of gas results in a gradual decline in both perigalacticon 
and apogalacticon and downwards jumps in apogalacticon near 
perigalacticon passages. Once the gas is stripped, the lack of in¬ 
teraction with the hot medium results in minimal change to the 
orbit, although once again the tidally disrupted dwarf lacks a 
clear centre, resulting in spurious jumps in perigalacticon and 
apogalacticon as the calculated centre rapidly shifts. 

These changes are above the accuracy expected to be deter¬ 
mined by Gaia surveys ( |Lux et al.|2010| approximately 14% ac¬ 
curate in determinations of peri- and apogalacticon), and hence 
the orbital changes due to ram pressure stripping is important in 
tracing the orbits of satellite galaxies and streams backwards in 
time. 
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[Fe/H] [Fe/H] 


(a) Pseudo-isothermal (b) NFW 

Fig. 12. [Fe/H] distribution of stars in the pseudo-isothermal simulations at the point of maximum stellar mass. Time of maximum stellar mass is 
shown. Stars formed before infall are shown in each plot in grey. 


The change in energy allows for the calculation of the coef¬ 
ficient of drag at each time step. As something potentially useful 
for semi-analytic models but not a focus of this paper, we dis¬ 
cuss this in Appendix]^ In general, the coefficient of drag, using 
a reference area of the region where Ahi >10 ^^ cm“^, is a factor 
of a few ~3-4 that of a sphere. 
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(a) Pseudo-isothermal (b) NFW 

Fig. 13. [Mg/Fe] of stars weighted by mass of each stellar particle in the pseudo-isothermal simulations at the point of maximum stellar mass. The 
black dashed line shows a running median over each 0.5 dex bin in [Fe/H]. 
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Fig. 14. Distribution of the radius of gas at infall inside the simulations, which goes on to form stars by the point of maximum stellar mass. 
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Fig. 15. Specific energy (kinetic plus potential energy) relative to infall and orbital evolution of the ram pressure models. RPS-50 is shown with 
a blue line, RPS-100 with red, NRPS-50 with cyan, NRPS-100 with magenta. For the orbital evolution, at each snapshot we calculate the peri- 
and apogalacticon of the orbit of the dwarf, based on its current position and energy. The galactocentric distance of the dwarf at any given time is 
shown with a solid line, while the peri- and apogalacticons are shown with dashed lines. 
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5. Conclusions 


For the first time, a full chemodynamical model of a dwarf 
galaxy is examined as it infalls into a the hot corona of another 
galaxy. Using an improved version of GEAR, we simulated the 
infall of a Sextans-like dwarf with a pseudo-isothermal and an 
NFW profile through utilisation of a rotating frame of reference 
with variable external gas density. The methods used here are 
novel and allow quick calculation of the ram pressure forces 
upon a dwarf. 


If the dwarf was already extended upon infall, e.g. hav¬ 
ing low-density cores due to early internal processes, even at 
large galactocentric distances the pressure of the hot halo more 
than doubles the star formation inside the dwarf. If the dwarf 
is denser, however, and retains an NFW profile, this effect is 
lessened although still noticeable. This star formation blows out 
large amounts of gas, greatly as sisting ram pressur e in stripping 
the dwarf, an effect also seen in |Gatto et al.| ( |2Q13] ). 

As gas is stripped from the dwarfs, it is pushed off centre 
against the direction of motion, resulting in a gradual energy 
loss of the dwarf altering the orbit. This gas is then finally lost 
around perigalacticon for most of the dwarfs considered here, 
after the dense medium causes a final burst of star formation. 
This final burst drags some dark matter with it due to the mutual 
gravitational attraction between the gas and dark matter. At peri¬ 
galacticon, this has the affect of greatly increasing the impact of 
tides, resulting in the pseudo-isothermal dwarf losing all its mass 
even at large perigalacticons (rperi~100 kpc), and transforming 
the dwarf into a stream. When gas is not lost at perigalacticon, it 
can be held off-centre until star formation eventually forces the 
gas outside of the dwarf, even as the ram pressure forces acting 
upon it drop. 

That the dense NFW models also succumb to these tidal ef¬ 
fects suggests that the Sextans dwarf spheroidal was markedly 
different from its present day form upon infall, and has been al¬ 
tered to its present day extended low-density form through the 
synergy of tidal and ram pressure forces. As the star formation 
is required to experience shut down within 6 Gyr, higher orbits 
are unlikely to solve these problems. 

The present day form of Sextans dSph could be achieved 
through it beginning with a denser dark matter profile at infall, or 
by having a much lower baryon fraction (well below the univer¬ 
sal ratio, e.g. through UV heating) to minimise these synergistic 
effects. 

The metallicity of the dwarf is also affected, with the rapid 
star formation being inefficient at recycling of metals, resulting 
in a lower peak in metallicity than would be expected by stellar 
mass alone. Although dwarfs are expected to be inefficient re¬ 
gardless ( |NRJ14| ), the evolution here is effectively halted for the 
pseudo-isothermal dwarfs. In addition to this poor recycling, the 
gas that goes on to form stars upon infall is not the same gas that 
would form stars in isolation. The combination of compression 
due to ram pressure and Bondi-Hoyle-like accretion results in a 
large fraction of gas from the outskirts going on to form stars. 
This star formation from metal poor gas makes dwarf galaxies 
whose kinematics suggest an early infall one of the best targets 
for searches of extremely metal poor stars. 

Even if the [Fe/H] evolution can continue, and in dwarfs 
that orbit at high-perigalacticon it may even be enhanced, the 
increase in star formation alters the ratio of Type II supernova to 
Type la supernova, resulting in slightly higher [Mg/Fe] at high 
[Fe/H]. This higher [Mg/Fe] is likely not directly observable, 
but notably, the dispersion is reduced with the ram pressure mod¬ 
els tending to be concentrated above [Mg/Fe] = 0, while the 
isolation models have a larger spread at high metallicity. 

When considering semi-analytic models of galaxy formation 
that does not include ram pressure, the most obvious changes to 
approximate ram pressure stripping are: 

- Dwarf galaxies that show signs of tidal disruption without 
considering ram pressure will likely not survive the combi¬ 
nation of tidal and ram pressure forces. 

- [Fe/H] evolution is slowed, moving a dwarf upwards on the 
[Fe/H] V L relation. 
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- At the high-metallicity end (past the “knee"), the [Mg/Fe] 
relation is marginally higher than would otherwise be ex¬ 
pected. 

- The orbital properties of the dwarf change, depending on 
the timescale of gas removal. Energy losses of 10-30% are 
readily possible, disproportionately resulting in a drop in 
apogalacticon. 

Slight errors in the ram pressure force are introduced 
throughout the calculation because of the method (10%-30% 
at perigalacticon, discussed further in appendix E- These er¬ 
rors are still well within the order of magnitude uncertainty in 
the observational values of the density and orbital velocity of 
the dwarfs at perigalacticon, and in particular seem to result in 
an underestimate of the ram pressure force at perigalacticon, as 
such any quantitative values (such as time of stripping) may be 
off. Nevertheless, the large impact of ram pressure stripping seen 
here suggests that these quantitative values are much better ap¬ 
proximations than those taken from dwarf galaxies modelled in 
isolation or in constant density wind tunnels and hence through 
the modelling of ram pressure stripped dwarfs, we can gain valu¬ 
able insight into the nearby satellite galaxies. 

References 

Agertz, O., et al. 2007, MNRAS, 380, 963 
Aoki, W., et al. 2009, A&A, 502, 569 

Arraki, K. S., Klypin, A., More, S., & Trujillo-Gomez, S. 2014, MNRAS, 438, 
1466 

Battaglia, G., Tolstoy, E., Helmi, A., Irwin, M., Parisi, R, Hill, V., & Jablonka, 
P.2011, MNRAS, 411, 1013 

Bland-Hawthorn, J., Sutherland, R., Agertz, O., & Moore, B. 2007, ApJ, 670, 
L109 

Boylan-Kolchin, M., Bullock, J. S., & Kaplinghat, M. 2011, MNRAS, 415, L40 
Breddels, M. A., & Helmi, A. 2013, A&A, 558, A35 

Chiboucas, K., Jacobs, B. A., Tully, R. B., & Karachentsev, I. D. 2013, AJ, 146, 
126 

de Blok, W. J. G., Walter, F., Brinks, E., Trachternach, C., Oh, S.-H., & Kenni- 
cutt, Jr., R. C. 2008, AJ, 136, 2648 

Gatto, A., Fraternali, R, Read, J. L, Marinacci, R, Lux, H., & Walch, S. 2013, 
MNRAS 

Governato, R, et al. 2010, Nature, 463, 203 
Grcevich, J., & Putman, M. E. 2009, ApJ, 696, 385 
Hopkins, P. R 2013, MNRAS, 428, 2840 
Huber, D., & Pfenniger, D. 2001, A&A, 374, 465 

Kafle, P. R., Sharma, S., Lewis, G. R, & Bland-Hawthorn, J. 2014, ApJ, 794, 59 
Kazantzidis, S., Lokas, E. L., & Mayer, L. 2013, ApJ, 764, L29 
Kirby, E. N., Cohen, J. G., Smith, G. H., Majewski, S. R., Sohn, S. T., & 
Guhathakurta, P. 2011, ApJ, 727, 79 

Klypin, A., Kravtsov, A. V., Valenzuela, O., & Prada, R 1999, ApJ, 522, 82 
Kroupa, P. 2001, MNRAS, 322, 231 
Lux, H., Read, J. L, & Lake, G. 2010, MNRAS, 406, 2312 
Maio, U., Dolag, K., Ciardi, B., & Tornatore, L. 2007, MNRAS, 379, 963 
Mayer, L., Mastropietro, C., Wadsley, J., Stadel, J., & Moore, B. 2006, MNRAS, 
369, 1021 

Mikhailov, M., & Freire, A. S. 2013, Powder Technology, 237, 432 
Miller, M. J., & Bregman, J. N. 2013, ApJ, 770, 118 

Moore, B., Ghigna, S., Governato, R, Lake, G., Quinn, T., Stadel, J., & Tozzi, P. 
1999, ApJ, 524, L19 

Navarro, J. R, Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 
Nichols, M., & Bland-Hawthorn, J. 2011, ApJ, 732, 17 
—. 2013, ApJ, 775, 97 

Nichols, M., Lin, D., & Bland-Hawthorn, J. 2012, ApJ, 748, 149 
Nichols, M., Revaz, Y., & Jablonka, P 2014, A&A, 564, A112 
Oh, S.-H., et al. 2015, AJ, 149, 180 

Ploeckinger, S., Hensler, G., Recchi, S., Mitchell, N., & Kroupa, P. 2014, MN¬ 
RAS, 437, 3980 

Pontzen, A., & Governato, R 2012, MNRAS, 421, 3464 

Revaz, Y, Arnaudon, A., Bonvin, V., Nichols, M., & Jablonka, P 2015, in prep 

Revaz, Y, & Jablonka, R 2012, A&A, 538, A82 

Rocha, M., Peter, A. H. G., & Bullock, J. 2012, MNRAS, 425, 231 

Roediger, E., et al. 2015, ApJ, 806, 104 

Sawala, T., Scannapieco, C., & White, S. 2012, MNRAS, 420, 1714 
Shetrone, M. D., Cote, P, & Sargent, W. L. W 2001, ApJ, 548, 592 

Article number, page 22 of[^ 


Smith, R., Fellhauer, M., & Assmann, P. 2012, MNRAS, 420, 1990 

Springel, V. 2005, MNRAS, 364, 1105 

Sutherland, R. S., & Dopita, M. A. 1993, ApJS, 88, 253 

Tafelmeyer, M., et al. 2010, A&A, 524, A58 

Theler, R., et al. 2015, in prep 

Torrey, P, et al. 2015, MNRAS, 447, 2753 

Walker, M. G., & Penarrubia, J. 2011, ApJ, 742, 20 

Wiersma, R. R C., Schaye, J., Theuns, T., Dalla Vecchia, C., & Tornatore, L. 

2009, MNRAS, 399, 574 
Xue, X. X., et al. 2008, ApJ, 684, 1143 
Zolotov, A., et al. 2012, ApJ, 761, 71 

Acknowledgements. The data reduction has been per¬ 
formed using the parallelised Python pNbody package 
(http://lastro.epfl.ch/projects/pNbody/). This work was sup¬ 
ported by the Swiss National Science Foundation (FNS 200021_153234). This 
work also benefitted from the International Space Science Institute (ISSI) in 
Bern, thanks to the funding of the team, “The first stars in dwarf galaxies". 


Appendix A: Rotation scheme integration and tests 

Because of the non-symplectic nature of the leapfrog integrator, 
care needs to be taken when dealing with the rotating reference 
frame to advance the velocity at each step due to the accelera¬ 
tions dependence on the velocity. 

We test a few schemes through two tests: first, a simple sys¬ 
tem of rotating rings, and also a dwarf that only experiences tidal 
forces as in lNRJ14] 

First, we consider a system of two simple rings, each com¬ 
prising a number of particles of negligible mass (we use 1000 
particles of mass 0.01 M©) in circular orbits at 50 kpc and 
100 kpc around a 1 x 10^ ^ M© point mass (modelled as a Plummer 
sphere with a softening of 0.5 kpc). If the central mass is con¬ 
sidered the origin, then the rotating box is positioned at (100,0) 
and is set to orbit in the same direction as the outer ring (clock¬ 
wise), giving the outer ring zero velocity in the rotating frame of 
reference. The inner ring is set to rotate in the opposite direction 
(anti-clockwise in this case). This system is then simulated for 
23.3 Gyr (5000 time units in our set up), just shy of 2.5 rotations 
of the outer ring. We show the results in Fig. |A.l[ We note that all 
errors here are time step dependent, with all methods perform¬ 
ing well at small time steps, however, large differences are seen 
at the time steps selected by the code. 

As the periodic box is in a circular orbit, this test only eval¬ 
uates the Coriolis and centrifugal forces. In the ideal case, the 
two rings (shown with black circles) maintains their orbits for all 
time, and this appears in the standard GEAR code (blue crosses; 
we note that this is identical to the standard Gadget-2 code, 
with an external potential force added) without the rotating ref¬ 
erence frame. Inside the rotating reference frame, calculating the 
forces at the standard time step (yellow squares) and using the 
half-time step velocity inside the fictitious forces gradually in¬ 
troduces errors that build in time, with the inner ring expanding 
to a radius of ~ 65 kpc and the outer ring staying constant (a sign 
that the centrifugal force is being evaluated accurately). We find 
better results (green points) if we calculate the fictitious forces 
using an approximated future velocity (Mid-L>) generated by the 
following: 

Vf,pred = yt-m + a(%_i/2, x^, V^_i/2, tU^_i/2)Ar/2, (A.l) 

V^+1/2 = V^_i/2 + a(q^_i/2, X^, v^,pred, u;^-i/ 2 )Al (A.2) 

Updating tu or q to an approximation does not greatly change 
the result here. In this case, the rings seem to suffer minimal 
contraction or expansion. 
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Fig. A.l. Results at 23.3 Gyr for the two-ring test. Two rings rotate in 
opposite directions around a 1 X 10^^ M© mass located at (0,0). When 
the frame of reference rotates (all but the standard GEAR code), it rotates 
around the point (100,0). The points should continue to orbit at 50 kpc 
and 100 kpc represented by two black rings. The standard GEAR code 
in a non-rotating frame is shown with blue crosses. In a rotating frame 
of reference, the results of calculating the forces at every time step is 
shown with yellow squares, using the half time step velocity (Mid-i>) 
as green points, and the iteration scheme (It-i;) with red triangles. Only 
every 20th point is shown for clarity. 

We display one further case (It-v, red triangles).This form 
performs better in other tests, although it is seemingly more sen¬ 
sitive to step size (at low-step size all methods perform similarly) 
and in pure dark matter simulations (which tend towards larger 
step sizes) should be avoided. Here, the mid point is approxi¬ 
mated as per above, but then it is approximated again using this 
calculated value. Here, we also assume that the angular velocity 
changes slightly across a time step, but that the acceleration in 


the angular velocity is negligible. We find 

Vr,pred = yt-1/2 + a(q^_i/ 2 , x^, V^_i/2, tu^_i/2)Ar/2, (A.3) 

^?,pred = + tfioldA^/2, (A.4) 

^f,pred ~ ^?-l/2 "f a(q^_l/29 X^, pj-ed? tL)^,pred)A^/2, (A.5) 

Vr+1/2 = ^t-l/2 + a(q^_i/ 2 , x^, V^,pred, OJt,^rQd)^t. (A.6) 


A second test is to compare the results with just tidal forces, 
particularly focusing on the star formati on prese nt within the 
dwarf. We follow the same method as in |NRJ14[ simulating a 
dwarf in isolation, before inserting it within the potential of a 
host galaxy. We choose an orbit with a perigalacticon of 80 kpc 
and an apogalacticon of 200 kpc, injecting the dwarf after 2 Gyr 
of evolution in isolation. This is the same simulated dwarf we 
use below, and its properties may be found in more detail in ^ 
The host halo parameters are the same as used in the main sim¬ 
ulations and are shown in Table The effects of just tides are 
discussed in detail in |NRJ14[ and here we just investigate the 
density profile, which provides resistance against ram pressure 
stripping, and star formation, which affects the chemical evolu¬ 
tion of the dwarf, to compare the codes. Here we implement the 



Fig. A.2. Stellar mass for the dwarf only experiencing tidal forces in the 
standard GEAR implementation (blue) and in the rotating reference with 
the Mid-i; scheme (green dash) and with the It-i; scheme (red dash-dot) 
and the It-t; scheme with a low step size (green dotted). The beginning of 
the orbit is shown with a vertical dotted line. For the first ~ 1. 5 Gyr , the 
stellar mass is basically the same with similar bursts (see Fig. |A.3| ), be¬ 
yond this, however, the Mid-i; scheme flattens out, still producing stars 
for another ~ 1.5 Gyr, but at a lower rate. The drop in stellar mass at 
the end in the rotation schemes is due to stellar particles being deleted 
as they cross the boundaries. This drop is lowe r in t he low step size due 
to better simulation of the dwarf core (see Fig. |A.4| ). 

particle deletion scheme discussed below in |2.3[ but this has 
minimal effects on the dynamics of the dwarf (we do not fol¬ 
low it all the way to destruction ). We co mpare the standard GEAR 
code (with slight updates from |NRJT^ , with the mid-point ve¬ 
locity scheme and the iterated velocity scheme for each property 
with the normal GEAR selected step sizes and the iterated veloc¬ 
ity scheme limited to the typical step size GEAR chooses for the 
main simulations (~0.094 Myr). 

We show the stellar mass in Fig. A.2| and the star forma¬ 
tion rate in Fig. |A.3| It is immediately obvious that the vari¬ 
ous schemes produce different star formation histories; while the 
standard GEAR code (top panels in both) experiences a period of 
extended star formation. In both of the other cases, this continues 
the bursty nature of the isolated dwarf (before the dashed line). 
That the GEAR code changes from a star formation history char¬ 
acterised by bursts to one with a gradual period of star formation 
(at least initially) is likely related to the lack of pause in cooling 
that occ urs when beginning a simulation from a snapshot (see 
NRJ14| for a discussion of this), why this does not occur with 
the other two is not known, but the star formation histories are 
similar apart from this initial region and that suggests it will not 
greatly affect the final results. 

A more noticeable change occurs in the density profile of 
the dwarf as shown in Fig. |A.4[ Here the mid-L> scheme (middle 
column) has a lower central density from the first perigalacticon 
passage (second row), which gets worse as time goes on. That 
much of the star formation occurs soon after input (and the dwarf 
is quenched in all cases) suggests that the schemes differ more 
in their estimation of the dynamical effects than in the chemi¬ 
cal evolution of the systems. The changes in the star formation 
(of order 20%) between the simulations is not a great concern 
when simulating dwarfs and is roughly the amount produced by 
changing random seeds. 
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Fig. A.4. Smoothed density in a central 1 kpc slice (logarithmic from maximum density at simulation start to 1% of this) through the centre of the 
dwarf in the three schemes. Each box is 10 kpc in width and height. The standard GEAR scheme is on the left, Mid-i> scheme in second column, 
It-v scheme with natural step sizes in the third, and the It-i; scheme with limited step sizes in the right column. From top to bottom, the slices take 
place at 2 Gyr (simulation start), 4.8 Gyr (perigalacticon), 7.6 Gyr (apogalacticon), and 10.4 Gyr (perigalacticon). In the rotating frame schemes, 
the dwarf is rotated back to the host galaxy frame of reference to be directly comparable with the standard GEAR scheme. As should be the case, the 
density profiles are identical at the simulation start, by the first perigalacticon, however (second row), the Mid-t; scheme is less dense in the centre 
and slightly over dense in the tidal tails; this central under density is carried over to the first apogalacticon (third row) before the Mid-i; dwarf is 
effectively destroyed at the next perigalacticon (bottom panel). The It-t; dwarf survives these passages in slightly better shape, particularly at the 
second perigalacticon where it still retains a slightly denser core. When limiting step sizes to those typical of the RPS simulations, the It-i; scheme 
performs quite well, possessing a similar core to the normal GEAR scheme at each point. 
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Appendix B: Particle deletion/creation scheme 



Fig. A.3. Star formation rate for the dwarf only experiencing tidal forces 
in the standard GEAR implementation (blue, top) and in the rotating ref¬ 
erence with the Mid-L> scheme (green, middle) and with the It-i> scheme 
(red, bottom). The beginning of the orbit is shown with a vertical dotted 
line. For the first ~ 1.5 Gyr, the bursts are similar between all cases with 
short bursts of ~ 4 X10""^ M© (and a slightly extended one in the normal 
GEAR code), which leads to a similar chemical evolution. Beyond this 
point however, Mid-f run drops to ~ 1 x 10“^ M© before becoming ex¬ 
tremely low levels around 6 Gyr. Suggesting that the change in density 
profile impacts the star formation rate, and subsequently the feedback 
and chemical evolution within the dwarf. 


The creation and deletion of particles is not without error. The 
creation of particles ignores things like pressure when deciding 
where particles are placed and their velocity, and the deletion 
relies on experimental values which are context dependent. 


The simplest test case to examine what errors we can expect 
is a simple wind tunnel. Even here however there are a number of 
free parameters to consider: the velocity, pressure, and internal 
energy of the particles all influence the flnal error. We consider 
a wind tunnel of length 50 (50 kpc in physical units) and inject 
particles with a velocity of 0.5 (103 km s“^). All particles have 
a density of 9 x 10“^ (3.7 x 10“^ cm“^) and an internal energy 
generated from a temperature of 2 x 10^ K. These conditions and 
velocity are similar to a dwarf orbiting at the outer edge of the 
halo. Although this box is not rotating, we delete particles if they 
come within 0.1 of the outer edge (this is smaller than the value 
used for the rotating box, but as discussed below, the rotation 
minimises the backwards wave). 

We displ ay th e velocity of the particles as a function of po¬ 
sition in Fig. Q and the ram pressure force divided by the true 
velocity in Fig. |B.2| at a time of 10. The creation and deletion 
of the particles introduce standing waves into the box, which are 
damped in amplitude towards the middle. Beyond the middle, 
the deletion of particles (and the subsequent vacuum it creates) 
increases the velocity beyond the centre. This creates a slight in¬ 
crease in velocity, which travels backwards before settling at a 
value of -0.55 + 0.02. If a deletion region of size 0.5 is used in 
the wind tunnel, the velocity instead settles at a value of -0.8. 
The change in velocity and changes in density due to the waves 
introduces a spread in ram pressure at the centre from about 90% 
of the true value to about 130%. 

We also consider the effects of the orbit changing, which in¬ 
troduces a pressure gradient across the box. The periodic box is 
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Fig. B.l. Two-dimensional histogram of the velocity and position of 
the particles in a wind tunnel (darker shades represent more particles 
there). The expected velocity is shown with a red dashed line at 0.5. 
The creation/deletion of particles introduces waves into the box, which 
are damped with position, and the deletion of the particles leads to a 
rising velocity towards the outer edge of the box. 



Fig. B.2. Two-dimensional histogram of the ram pressure force (di¬ 
vided by the expected) and position of the particles in a wind tunnel 
(darker shades represent more particles there). The expected ram pres¬ 
sure force is shown with a red dashed line at 1.0. Here, we see that 
the creation/deletion of particles introduces waves into the box, with a 
rising velocity towards the outer edge of the box. 


placed in an orbit with a perigalacticon of 80 kpc and apogalacti- 
con of 200 kpc, identical to that used in Appendix [A| in a box of 
side 50 kpc, a gas particle mass of 6.0 x 10^ Mq, and a halo den¬ 
sity profile as per |3.2| We show the results in Fig. |B.3| as a func¬ 
tion of the radius of the orbit. Here, the pressure gradient results 
in a slight underestimation of the ram pressure at perigalacticon 
(by about -10%), while being slightly above at apogalacticon. It 
is clear that over the densest fastest region of the orbit, any dwarf 
likely experiences a slight underestimation of the ram pressure, 
with the overestimation that occurs towards apogalacticon. 

Notably, the larger error appears to be as the dwarf increases 
in galactocentric distance (i.e. approaching apogalacticon), and 
suggests that the scheme has more trouble accounting for a de¬ 
crease in pressure than the increase in pressure associated with 
infall. 
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Fig. B.3. Ratio of the ram pressure force near the centre of the box to 
the expected versus radius for a sample orbit. The region of the orbit 
decreasing in galactocentric distance is shown in blue, while the region 
of increasing distance is shown in red. 

The choice to base it on the current positions of the dwarf 
results in creating an effective delay in the density of the halo 
reaching the dwarf (since the particles have to travel half the 
box). As the orbit can change on these scales, the solution to this 
problem would require integrating the pivot point forward and 
inserting the particles at its future density. Given the correction 
factor, this is not a trivial solution, particularly at perigalacticon 
where the impact is biggest, however, future versions of the code 
will attempt to incorporate this solution. 

Appendix C: Particle creation 

At simulation start, the positions of particles of an SPH glass (of 
A^giass particles) in a unit cube are passed through to memory. At 
each time step, of size At, particles of mass mgas are injected with 
a density corresponding to the halo halo at density phaio- Particles 
are inserted by mapping the volume required to be filled by the 
particles to the glass cube. We show a schematic of the particle 
creation in Fig. |C.l| to assist in the understanding. 

If the number of particles required to fill the entire box (of 
side L) is Ahaio = Phaio^^^gas. then the volume at each time step 
that needs to be filled corresponds to a rectangular prism (xxyxz) 
on the glass of size, 

A^halo IIVpIlAf / jVhalo f / jVhalo 

glass / L \N glass / \^glass / 

All particles within this prism, and moving the prism along the 
velocity axis at each time step, are “accreted", that is, a new 
gas particle is created at the position (A/haio/A^giass)~^^^(-^^, yg, Zg), 
where (Xg, ijg, Zg) is the position of the particle on the glass rel¬ 
ative to the origin of the prism. As (A^haio/^^giass)”^^^ is not, in 
general, equal to one the periodicity is lost, but we find using a 
glass is still considerably better than randomly distributing halo 
particles in the new volume. 
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Periodic Box 



Glass 



Fig. C.l. Particle creation schematic. The periodic box is shown in the left panel. The SPH glass is shown on the right. When particles need to 
be created (below the dashed line), a region of the SPH glass is selected (right panel, grey shading) according to equation ( |C.1| ). The positions of 
the particles within this region is then scaled to match the boundaries of the particle creation region, resulting in the grey points. The velocity of 
the particles is then set according to equation ( p^ , and the other properties based on the assumed halo profile. At the next time step, the region of 
glass selected would start at the top of the shaded region (0.1 in the right-hand box) repeating the glass in any dimension when necessary. 


Appendix D: Comparison of baryon free sateiiites 
and haios 

We compare here the mass and orbital evolution of a satellite 
with dark matter only and one with baryons and dark matter 
that orbits in the rotating reference frame without the external 
host coronae being present. The baryon free model, NBF-50, 
is created by evolving the same original dark matter potential 
without dwarfs and inserting it on the same orbit as NRPS-50, 
which has a perigalacticon of 50 kpc and an apogalacticon of 
200 kpc. The tidal model, NTO-50, is simply the NISO model 
inserted at 2 Gyr on the same orbit of 50 kpc perigalacticon and 
200 kpc apogalacticon, but with no Milky Way-like hot halo. In 
both cases, the models were limited to a maximum step size of 
~0.094 Myr, and the step size GEAR quickly selects for the ram 
pressure simulations; all other parameters are the same as for the 
RPS cases. 

These two models provide complementary tests. NBF-50 en¬ 
sures any orbital evolution is not due to dynamical friction or 
errors in the particle creation/deletion, creating spurious gravita¬ 
tional attractions behind the dwarf. NTO-50 checks that the mass 
evolution is not a by-product of the rotation scheme or only due 
to tidal influences. Both models, without any form of hydrody¬ 
namic drag, also provide a test that the orbital evolution is not a 
by-product of the rotation scheme. 

Appendix D. 1: Mass evolution 

We show the dark matter mass of the NTO-50 and NBF-50 mod- 
els and NRPS-5Q as a function of time in Fig. |D.1| Similar to 
Arraki et al.| ( [2Q14| ), the removal of baryons in NBF-50 pro¬ 
duces a dramatic loss in dark matter density (and therefore mass 
within the 1 kpc and 3 kpc boundaries), although the mass loss is 


slightly under that found by [Arraki et al. [ possibly attributable to 
the instantaneous removal of the baryons here. The tidal model 
NTO-50 retains the central mass density of NRPS-50, and tracks 
it closely until the first perigalacticon. 

At the first perigalacticon, mass is lost from ah models, how¬ 
ever, it is dramatically more from NRPS-50 than NBF-50 or 
NTO-50, which lose ~30%-50% within 3 kpc. This is com¬ 
pared to NRPS-50 which because of the ram pressure loses more 
than 70%. By the next perigalacticon, the NRPS-50 model is de¬ 
stroyed into a stream, while both the NBF-50 and NTO-50 mod¬ 
els survive until the present day. 

That both models undergo tidal stripping but retain their gen¬ 
eral structure confirms that the ram pressure stripping of gas 
acts synergisticahy with tidal stripping to remove mass from the 
dwarf, including its dark matter mass. 

Appendix D.2: Orbital evolution 

We compare t he orb ital evolution (see |4.6[ ) between the three 
models in Fig. [D.2[ Notably, neither NBF-50 or NTO-50 show 
a change in orbital energy, as these orbital energies are always 
approximately the equal to the infah energy. Similarly, their 
peri- and apogalacticons remain approximately constant. NBF- 
50 does show a minor drop in apogalacticon of order 5 kpc, 
which can be regarded as the error arising because of the halo 
generation method. Notably, this is still much smaller than the 
50 kpc drop in apogalacticon NRPS-50 suffers at the first peri¬ 
galacticon. The constant nature of these peri- and apogalacticons 
increases the confidence that the drag experienced on the dwarf 
is due to a hydrodynamic interaction between the gas of a dwarf 
and the hot halo, which is then gravitationally coupled to the 
dark matter of the dwarf. 


Article number, page 27 of[^ 












A&A proofs: manuscript no. rps 



Time [Gyr] 


Fig. D.l. Mass of dark matter within 1 kpc (thin lines) and 3 kpc (thick 
lines) for the NRPS-50 (cyan), NBF-50 (brown), and NTO-50 (pink) 
models. 



Fig. D.2. Orbital evolution of NRPS-50 (cyan), NBF-50 (brown), and 
NTO-50 (pink). Line styles are as per Fig.p~5] 


Appendix E: Drag of a dwarf galaxy 

The drop in energy throughout a dwarf’s orbit (Fig. [Tg and sub¬ 
sequent orbital change suggests that hydrodynamic drag is non- 
negligible when considering the evolution of a dwarf. By com¬ 
paring the energy change to the mass nearby the coefficient of 
drag of the galaxy can be calculated. If the kinetic energy at time 
step i is and potential energy then the coefficient of drag 
at the time step i is given by 

Cdi = — - ^^(E.l) 

PiAivfMi(ti+i - li) 

We consider the mass of the dwarf within 10 kpc of the cen¬ 
tre, and assume the energy change is simply the energy per unit 
mass of the centre multiplied by the mass. The choice of area 
is much more ambiguous, consisting of a clumpy medium, and 
it is not clear what reference should be chosen as hot gas may 
stream through parts of the structure. At each snapshot we select 
the total area of the dwarf, which exceeds a column density in 
the direction of motion of = 1 x 10^^ cm“^ averaged over 
0.25 kpc^ cells. This boundary is not the only potential choice, 
but seems to correspond to the area of the dwarf through which 
hot gas cannot easily flow. 

As the coefficient of drag generally depends on the Reynolds 
number of the flow, we calculate the Reynolds numbeij^at each 
time step as well, assuming the reference length is the same as 
the radius of a circle with equivalent area to the reference area. 
We plot the mean coefficient of drag versus Reynolds number in 
Fig. |E.l| for the NRPS-50 model for the first 1.5 Gyr (450 snap¬ 
shots) and the NRPS-100 model for ~8 Gyr (2000 snapshots). 
This time selection avoids the recreation event in NRPS-50, 
which leads to significant errors and stream formation in NRPS- 
100. We also show a fit to the coefficient of drag to a smooth 
sphere versu s Reynolds number (equation 7 of [Mikhailov &| 
|Freire|20T3] ) and see that although the coefficient of drag on the 
dwarf is much higher (as could be expected from the shape of the 
dwarf), it follows a similar trait of rapidly dropping until Re~200 
where Cd~3-10 before dropping down to Cd~l. 

The large amount of scatter at low Reynolds number is not in 
itself surprising, with the dwarf changing shape numerous times 
due to internal feedback, which is minimised at high Reynolds 
number because of the higher containing pressure. This large de¬ 
gree of scatter (and the dependence of Cd on area) makes it diffi¬ 
cult to approximate the orbital change in semi-analytic models. 


^ We note that the artificial viscosity affects the Reynolds number of 
the flow and is different between each particle pair. We ignore these 
differences and use the viscosity of the hot gas halo rather than the nu¬ 
merical viscosity. 
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Fig. E.I. Coefficient of drag of a dwarf galaxy throughout its orbit 
against the Reynolds number of the flow (calculated using the hot gas 
viscosity, not the numerical). The solid lines show the mean of the co¬ 
efficient of drag for NRPS-50 (cyan) and NRPS-100 (magenta) within 
each (logarithmic) Reynolds number bin. Dashed lines show the stan¬ 
dard deviation on each bin. A fit to the changing coefficient of drag of 
a sphere versus Reynolds number from [Mikhailov & Freire| ( |2013j ) is 
shown with a black line. 
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